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Figure  1:  Mesa,  Shelf  and 
Pedestal  eonstraints  are  chosen 
to  optimize  an  array  factor  with 
this  sombrero  shape. 


1  Introduction 

Suppose  we  have  a  planar  transmit  array  of  weighted  antenna  elements  laid  out  on  a  lattice, 
a  regular  grid,  and  suppose  each  embedded  element  pattern  is  identical  (typically  achieved 
using  guard  elements).  Then  the  far-held  array  pattern  is  just  the  product  of  the  embedded 
element  pattern  and  an  array  factor  that  is  the  2D  Fourier  transform  of  the  array  taper. 
Mathematically,  the  array  taper  comprises  impulses  at  the  element  locations  with  areas 
given  by  the  element  weights.  We  can  speak  of  designing  the  weights,  the  taper,  or  the 
array  factor.  These  are  equivalent. 

Is  it  possible  to  optimize  a  transmit  array  taper  that  illuminates  both  a  highly  directive 
cluster  of  receive  beams  and  a  broad  region  of  the  sky,  perhaps  to  support  a  ubiquitous  radar 
function  E 13?  Such  a  transmit  array  factor  might  allow  us  to  track  several  targets  and 
simultaneously  monitor  a  broad  region  of  interest  for  new  targets.  The  goal  of  this  report  is 
to  determine  whether  it  is  feasible  to  design  a  transmit  array  taper  with  such  characteristics. 
The  larger  issues  involved  in  designing  such  a  radar  system  are  outside  the  scope  of  this 
report  and  will  not  be  addressed. 

The  array  factor  that  we  desire  has  what  we  call  a  sombrero  structure  (Figure [T]),  where 
the  term  mesa  refers  to  the  top  of  the  sombrero,  shelf  xtftxs.  to  the  top  of  the  sombrero’s 
brim  and  pedestal  refers  to  the  bottom  of  the  sombrero’s  brim.  The  optimized  amplitude- 
taper  designs  that  follow  have  constraints  in  one  or  more  of  the  angle  regions  defined  by 
these  terms.  Constraints  on  the  mesa  region  limit  passband  ripple,  constraints  on  the  shelf 
region  limit  the  array-factor  magnitude  in  the  stopband,  and  constraints  on  the  pedestal 
region  define  the  height  of  the  sombrero’s  brim. 

In  this  report  we  are  particularly  interested  in  an  “active”  transmit  array,  in  which  each 
element  is  driven  by  its  own  power  amplifier.  There  are  two  basic  categories.  In  one  cat¬ 
egory,  the  power  radiated  is  limited  by  the  power  available  from  individual  elements.  The 
output-power  capabilities  of  the  individual  amplifiers  driving  the  elements  is  the  determin¬ 
ing  factor.  In  the  other  category,  the  power  radiated  is  limited  only  by  a  constraint  on  total 
transmit  power  as  determined  by  available  prime  power  or  by  heat-dissipation  issues,  and 
individual  amplifiers  can  provide  power  as  needed.  These  two  categories  are  idealizations; 
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in  real  systems  power  might  well  be  constrained  by  both  a  total-power  limit  and  by  per- 
element  power.  However,  consideration  of  these  two  extremes  reveals  a  key  choice:  given 
a  desired  target  illumination  level,  should  our  weight  design  minimize  peak  element  power 
or  total  transmit  power? 

Minimizing  peak  element  power  (without  other  design  constraints)  leads  naturally  to 
a  phase-only  taper.  All  element  weights  have  the  same  magnitude.  To  minimize  total 
transmit  power,  however,  complex  weights  are  needed,  and  design  symmetries  typically 
lead  to  weight  phases  being  used  only  for  beam  steering  so  that  the  unsteered  taper  is  a 
pure  amplitude  taper. 

Amplitude  tapers  are  our  main  focus,  not  because  they  are  what  is  most  needed  or 
desired  but  simply  because  we  have  first-rate  tools  for  designing  complex  or  amplitude 
tapers  but  are  limited  to  the  crudest  of  design  methods  for  phase-only  tapers.  That  situation 
is  not  entirely  about  software  and  methods  but  is  about  the  problem  itself:  a  phase  is  one 
degree  of  freedom,  but  a  complex  weight  is  two.  Should  we  not  expect  that  doubling  the 
number  of  degrees  of  freedom  will  permit  better  optimization?  Further  and  more  important 
is  that  the  two  design  spaces  are  different  in  a  fundamental  geometric  sense.  The  amplitude- 
taper  space  is  smoothly  shaped  and  easily  searched,  while  the  phase-only  space  is  jagged 
and  rough,  a  problem  only  made  incomprehensibly  worse  by  the  high  dimensionality — 
often  thousands — of  the  space.  It’s  largely  unsearchable  except  by  trial  and  error. 

So  here  we  explore  complex  or  amplitude  tapers  at  some  length  and  phase-only  tapers 
quite  tentatively.  Phase-only  tapers  are  a  subset  of  complex  tapers  of  course,  so  if  there  is 
no  good  solution  to  the  complex-taper  sombrero  design  problem,  there  cannot  be  a  good 
solution  to  the  phase-only  sombrero  problem  either.  Addressing  complex  tapers  first  then 
implicitly  seeks  to  learn  whether  there  is  something  in  the  nature  of  array  factors  themselves 
that  precludes  a  sombrero  structure.  If  our  answer  had  been  ”no  sombreros  allowed”,  it 
would  have  applied  to  both  types  of  tapers.  As  it  turns  out,  that  is  not  our  answer. 

Certain  details  are  common  to  every  design  presented  in  this  report.  All  of  the  designs 
in  this  report  assume  a  roughly  circular  4507  element  array  arranged  on  a  A/v^  spaced 
equilateral  triangular  grid  (Figure  |^.  For  ease  of  optimization,  we  require  that  the  taper 
be  unchanged  when  the  whole  array  is  subjected  to  any  of  the  five  indicated  rotations  or 
mirrored  through  any  of  the  six  lines  shown  in  Figure]^  This  automatically  makes  the  array 
factor  invariant  to  those  same  rotations  and  reflections.  These  array  symmetries  reduce  the 
number  of  optimization  variables  by  a  factor  of  roughly  12.  A  real  design  might  well 
call  for  a  nonsymmetric  pattern,  but  this  approach  provides  the  essential  insight  into  the 
trade  space  at  a  much  reduced  computational  cost.  We  assume  a  narrowband  array  and  do 
not  consider  the  wideband  case.  In  each  design  we  show  boresight  beams,  however  it  is 
possible  to  phase  steer  any  of  these  array  factors  off  of  boresight. 

Section  1^  presents  optimized  amplitude  taper  designs.  These  amplitude  taper  designs 
start  with  a  very  simple  pencil  beam  design  and  work  towards  the  sombrero  pattern  de¬ 
scribed  above.  The  report  includes  a  progression  of  designs  to  show  the  effects  of  each 
constraint  that  is  added  to  the  optimization  to  eventually  produce  the  sombrero  pattern. 

Section  [^presents  several  phase-only  tapers  constructed  as  2D  chirps,  with  linear-chirp 
designs  followed  by  nonlinear-chirp  designs.  Chirp  parameters  were  chosen  and  tweaked 
by  hand,  simply  by  trial  and  error.  The  results  are  rough  and  serve  mostly  to  illustrate  the 
severity  of  the  phase-only  limitation. 
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Figure  2:  Element  layout 


Figure  3:  The  optimized  tapers 
(and  array  factors)  are  made 
invariant  to  various  reflections 
and  rotations  simply  to  make 
optimizing  a  large  array  man¬ 
ageable. 


2  Amplitude  Tapers 

In  a  well-matched  transmit  array  the  power  into  each  element  is  proportional  to  the  squared 
magnitude  of  the  corresponding  element  weight,  so  the  total  power  input  to  the  array  is 
proportional  to  the  sum  of  the  squared  weight  magnitudes,  the  energy  in  the  weight  func¬ 
tion.  In  array  optimization  we  are  generally  interested  in  maximizing  the  power  density 
at  some  point  in  the  far  field,  but  should  we  do  this  assuming  that  the  total  power  deliv¬ 
ered  to  the  elements  is  fixed  or  that  there  is  some  limit  to  the  power  that  can  be  applied  to 
any  individual  element?  The  answer  depends  on  the  hardware  design,  and  so  each  design 
was  optimized  both  ways.  For  reasons  of  convenience  we  actually  fix  the  far-field  power 
density  at  some  point  in  space  and  minimize  one  of  these  two  measures  of  element  power. 
Fixing  the  far-field  power  density  and  minimizing  total/peak-element  power  is  equivalent, 
to  within  a  common  scale  factor  on  the  weights,  to  fixing  the  total/peak-element  power  and 
maximizing  the  far-field  power  density. 

We  use  our  MATLAB  toolbox  Opt  [[3]]  to  design  the  amplitude  tapers  that  follow.  The 
Opt  toolbox  solves  a  particular  class  of  convex-optimization  problems  in  which  a  linear 
objective  function  is  minimized  subject  to  any  number  of  constraints,  each  either  a  linear 
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constraint  or  a  convex  quadratic  constraint. 

The  first  two  designs  in  this  section  are  reference  cases,  the  Uniform  design  (page 
and  the  Mesa  Reference  design  (page[^.  Following  these  reference  cases  are  two  series  of 
optimized  designs:  designs  minimizing  the  energy  in  the  element- weight  function  ([Design 
|2.3|on  page[T^through|Design  2.9|on  page|^  and  designs  minimizing  the  peak  magnitude 
of  the  element- weight  function  ([Design  2.10|on  page  [T7[  through  [Design  2.16[on  page  [39]). 
For  each  series  of  amplitude  taper  optimizations,  we  include 

•  a  simple,  Taylor-like  pencil  beam  design, 

•  a  design  with  mesa  constraints, 

•  a  design  combining  the  mesa  constraints  with  a  Taylor-like  shelf  constraint, 

•  a  design  combining  the  Taylor-like  pencil  beam  and  the  brim  constraints,  and 

•  our  final  sombrero  designs. 


Design  2.1  A  simple  reference  case:  test  the  code 

The  reference  design  features  uniform  illumination.  Each  of  the  N  weights  is  set  to  1/A  so 
that  the  array  factor  is  unity,  0  dB,  at  boresight.  No  optimization  is  involved,  but  the  code 
for  computing  and  plotting  the  array  factors  and  performance  numbers  is  tested  in  a  simple 
way  here. 

This  reference  design  serves  a  second  purpose  as  well,  one  based  on  it  possessing  a  very 
special  property.  This  reference  design  can  be  represented  by  2D  discrete-time  Fourier  pair 
Tn  ^  ^(f).  where  taper  Vn  is  a  function  of  an  integer  two- vector  element  index  n.  In 
particular,  Vn  =  1/ A  for  n  in  the  array  and  zero  elsewhere.  Using  this  reference  design  the 
signal  power  radiated  in  the  boresight  direction  is  proportional  to  | i?(0)  p  =  1,  and  the  total 
power  delivered  to  elements  by  the  transmitter  is  proportional  to  ||r|p  =  knP  =  1/A. 
Classically,  boresight  gain  is  taken  as  the  ratio  of  these  quantities  |i?(0)p/||r|p  =  A. 
Likewise,  if  taper  hn  ^  H{i)  is  used  instead,  its  gain  relative  to  that  reference  gain  is 


|g(0)l7llftlP  „ 

|ii(0)|V||r|P  ||ft||2]V 


n 


||A|P  N 


\{h,r)\^ 

ll^lPIkP 


(1) 


using  on  the  right  inner  product  {h,  r)  =  h^r^.  The  point  of  this  seemingly  odd 

formulation  of  the  gain  improvement  is  that  the  ratio  on  the  right  cannot  exceed  unity, 
by  the  famous  inequality  of  Cauchy,  Schwarz,  and  Bunyakovsky.  Therefore,  no  taper  can 
improve  on  the  boresight  gain  of  this  reference  case,  so  for  each  design  to  follow  we  give 
the  dB  distance  below  this  gain  bound,  classically  referred  to  as  taper  loss,  as 


(weight-energy  taper  loss)  =  —10  login 


|(/i,r)|2 


(2) 


This  reference  case  also  serves  as  the  gain  bound  in  the  peak-element  power  sense.  If 
we  replace  the  element  energy  ||^|p  with  the  peak  element  magnitude  |/imaxp, 


|i^(0)|V|/^maxP 

I^WIVkmaxP 


n 


(3) 


But  the  numerator  is  equal  to  (/i, r)^/|rmaxp.  so  for  each  design  we  can  report  the  dB 
distance  below  the  gain  bound  as 


(max-weight  taper  loss) 


-lOlogio 


\{h,r)\^ 

|/imaxPiV2|r„,axP 


(4) 


Array  factors  are  plotted  in  what  is  often  and  somewhat  confusingly  called  either  sine 
space  or  direction-cosine  space.  The  far-field  plane  wave  radiated  with  3D  radian  vector 
wavenumber  27rk  has  amplitude  proportional  to  $(k)i/((kA  —  s)^B),  where  $(k)  is  the 
embedded  element  pattern  and  H  ((kA  —  s)^B)  is  the  phase  steered  array  factor  of  interest. 
Here  H{i)  is  the  2D  discrete-time  Fourier  transform  of  the  element  weights  hn,  s  is  a  unit 
vector  in  the  steering  direction,  and  B  is  a  dimensionless  3x2  matrix  characterizing  the 
geometry  of  the  element-layout  grid.  For  a  plane  wave,  normalized  wavenumber  kA  is  a 
unit  vector,  but  difference  kA  —  s  is  not  and  can  have  a  vector  length  approaching  two.  For 
this  reason,  we  plot  i/(k^AB)  without  restriction  to  the  unit-vector  case.  Decomposition 
k  =  k  +  kx  into  array  plane  and  perpendicular  boresight  components  allows  us  to  take 
advantage  of  the  orthogonality  of  k^  to  the  columns  of  B,  both  of  which  always  lie  in 
the  array  plane,  through  the  resulting  fact  that  k^B  =  k^B.  To  plot  201ogio 
then,  we  will  plot  20  log^g  |  versus  components  of  k  in  the  azimuth  and  elevation 

directions,  and  most  often  we  will  write  the  array  factor  as  |i7(k^B)  |  as  well. 

The  array  factor  is  displayed  at  three  different  zoom  levels  in  the  first  column  of  each 
figure. 


galaxy  zoom:  A  view  of  the  array  factor  showing  several  periods. 


world  zoom:  A  big-picture  view  of  the  array  factor.  The  latitude  and  longitude  lines  repre¬ 
sent  elevation  and  azimuth  respectively,  and  the  interior  of  the  “world”  is  the  visible 
region,  the  part  of  the  array  factor  corresponding  to  physical  directions.  Elevation 
and  azimuth  lines  are  spaced  at  15°  intervals. 


nation  zoom  60  dB:  A  close-up  view  of  the  array  factor  showing  details  of  the  transition 
region.  The  plot  spans  —7.5°  to  7.5°  in  elevation  and  azimuth,  with  a  color  scale  of 
60  dB.  Elevation  and  azimuth  lines  are  spaced  at  1°  intervals.  The  dB  color  scale  uses 
contrasting  colors  for  the  —  3  dB  and  —  2  dB  levels  in  order  to  make  the  associated 
beamwidths  visible. 


nation  zoom  4  dB:  A  close-up  view  of  the  array  factor  showing  details  of  the  mesa,  with 
a  color  scale  of  5  dB.  The  dB  color  scale  uses  contrasting  colors  for  the  —  3  dB  and 
— 2  dB  levels  in  order  to  make  the  associated  beamwidths  visible. 
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slice  plot:  An  array-factor  slice  plotted  on  both  a  linear  and  dB  scale.  The  slice  shown  is 
from  boresight  to  the  edge  of  the  visible  region  along  elevation.  This  plot  depicts 
qualitative  behavior  of  the  array  factor  and  is  particularly  helpful  when  the  optimiza¬ 
tion  behavior  is  unexpected. 

element  weight  vs.  radius:  The  element  weights  are  plotted  versus  the  radius.  Due  to 
the  symmetry  of  the  problem  specifications,  none  of  the  optimal  amplitude  weights 
have  nonzero  imaginary  components,  so  all  phases  are  0°  or  180°.  Negative  and 
nonnegative  element  weights  are  different  colors  (or  will  be  for  designs  with  weights 
of  both  polarities). 

element  weights:  The  element  weights  are  plotted  in  3D  perspective.  (This  plot  is  not 
particularly  revealing  for  this  design,  but  will  be  for  others  to  follow.) 

The  reported  taper  loss  and  taper  efficiency  numbers  are  two  different  ways  of  expressing 
the  same  measurement.  They  are  not  two  different  measurements. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency)  0  dB  (100%) 
max-weight  taper  loss  (max-weight  taper  efficiency)  0  dB  (100%) 


Figure  4:  Uniform  element  weights 
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Design  2.2  A  mesa  reference  case 

The  taper  losses  of  @  and  @  jump  upward  when  we  add  a  mesa  to  the  designs  because  of 
the  large  increase  in  noise-gain  factor 


E 


\K\^=  [  |i/(f)prff  =  A2|B^B|^/2  /  \H{k^XB)\'^dk, 

J  unit  J  one 


'unit 

square 


one 
period 


where  di  is  differential  area  in  two- vector  f  and  dk  is  differential  area  in  the  array-plane 
component  k  of  three-vector  k.  The  integral  on  the  right  becomes  much  larger  when 
H(^XB)  is  given  a  near- unity  value  over  the  entire  mesa  rather  than  over  a  pencil  beam 
as  before.  What  this  means  is  that  the  [Design  2.1  [reference  is  a  poor  choice  for  tapers  with 
mesas:  How  do  you  know  if  a  taper  loss  of  12  dB  is  good  or  not?  It  turns  out  that  we  can 
extend  the  notion  of  taper  loss  to  explicitly  account  for  a  mesa  in  a  way  that  is  particularly 
clean  for  amplitude  tapers. 

How  do  we  characterize  the  mesa  mathematically?  Array  factors  are  plotted  as  func¬ 
tions  of  the  array-plane  component  k  of  three-vector  k,  so  the  mesa  is  actually  a  set  of 
array-plane  values.  Each  such  value  can  be  written  as  a  weighted  combination  of  a  pair  of 
linearly  independent  array-plane  basis  vectors.  The  columns  of  basis  matrix  B  are  such  a 
pair,  but  it  is  convenient  here  to  instead  use  the  rows  of  B"^,  the  Moore-Penrose  pseudoin¬ 
verse  of  B  that  is  formally  given  as  B"*"  =  (B^B)“^B^and  that  is  realized  in  matlab  with  the 
pinv  ( )  function.  Its  rows  are  in  the  array  plane  because  each  is  just  a  linear  combination 
of  the  basis-vector  columns  of  B,  and  this  new  basis  is  appealing  because  if  we  represent 
an  array-plane  value  k  using  a  row  two-vector  f  of  coordinates  so  that  k^  =  yfB"^,  the 
simplifying  relationship  k^AB  =  fB'''B  =  f  results. 

Choosing  the  mesa  now  amounts  to  choosing  a  set  of  f  vectors.  We  will  soon  be  tak¬ 
ing  averages  across  the  mesa,  and  averages  are  nicely  handled  by  the  expected  values  of 
probability  theory,  so  let  us  suppose  that  f ,  and  implicitly  k  as  well,  is  a  random  vector  dis¬ 
tributed  probabilistically  across  one  array-factor  period  and,  in  particular,  with  probability 
concentrated  in  the  mesa.  Then 


average 
amplitude  gain 


E|H-(k'^AB)|  E\H{i)\ 


By  distributing  probability  in  the  mesa  according  to  the  weighting  desired  in  the  average, 
we  could  in  principle  obtain  a  weighted  average.  We  do  not  do  so,  however. 

Consider  some  specific  cases. 


Case:  single-point  mesa.  If  f  =  0  with  probability  one. 


average 
amplitude  gain 


|g(0)l 

ll^ll 


For  this  degenerate  mesa  the  average  amplitude  gain  is  just  the  Foresight  gain  used  to 
derive  weight-energy  taper  loss  in  [Design  2.1  [(which  begins  on  pageQ.  There  is  no 
reason  to  compute  average  amplitude  gain  separately  with  the  single  point  considered 
as  a  mesa. 
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Case:  mesa  defined  by  a  pdf.  This  is  the  main  case  of  interest.  Suppose  the  probability 
of  f  falling  into  a  subset  T  of  the  unit  square  is 

P[fG.F]  =^p(f)c?f, 

where  di  is  differential  area  in  the  f  plane  so  that  p(f)  on  the  unit  square  is  a  proba¬ 
bility  density  function  (pdf)  for  f.  Then  let  p(f)  be  periodic  with  the  unit  square  as 
its  period,  so  that  it  can  be  inverse  Fourier  transformed.  After  inverse  transforming 
|P'(f)  I  as  well,  we  can  use  2D  Fourier  pairs  hn  ^  \H{i)  \  and  p(f )  along  with 
Parseval’s  relation  and  ||(i||^  =  ||^||^,  which  is  obvious  in  the  frequency  domain,  to 
write 

arnpmudfgain  =  ^  (f ) |  p(f )  di  (5) 

square 

m  ■ 


There  are  two  subcases. 

Subcase:  zero-phase  taper.  The  amplitude  tapers  of  this  report  not  only  give  each 
weight  in  the  taper  a  zero  phase,  the  design  symmetries  imposed  also  give  each 
point  of  the  associated  array  factor  a  zero  phase.  Though  somewhat  confusing, 
the  custom  in  signal  processing  is  to  refer  to  tapers  with  the  latter  property  as 
zero-phase  tapers  (or  filters)  or,  less  precisely,  as  linear-phase  tapers  (or  filters). 
Typically  a  zero-phase  taper  not  only  has  H{i)  real  everywhere  in  the  mesa, 
generally  in  the  mesa  it  also  has  i7(f)  >  0  so  that  (f) |  =  ^(f)  and 

average  ^  1  /  u(c\  pc  ^ 

amplitude  gain  “  1 1  1 1  1  ^ 1 1  . 

square 

Since  fin  =  0  for  n  not  corresponding  to  actual  array  elements,  inner  product 
(fi,  (p)  =  fin^n  has  many  zero  terms.  We  can  just  as  well  decompose  4>  as 
0n  =  +  4’±n,  whcrc  and  can  be  nonzero  only  for  n  corresponding  to 

elements  actually  inside  and  outside  the  array  respectively,  and  write  (fi,  (p)  = 
(fi,  ^).  It  is  then  useful  to  consider  the  ratio 

average 

amplitude  gain  _  (fi, 

m 

because  by  the  Cauchy-Schwarz-Bunyakovsky  inequality  the  quantity  on  the 
right  will  never  exceed  one,  but  we  can  make  it  equal  to  one  by  setting  fi  =  ^. 
We  characterize  the  closeness  of  a  mesa  design  to  this  optimum,  using  termi¬ 
nology  loosely  (because  technically  we  are  averaging  the  gain  and  not  the  loss). 
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by  reporting 


average 


(average  mesa  taper  loss) 


(6) 


with  confidence  that  in  this  zero-phase  subcase  it  must  be  nonnegative. 

Subcase:  a  taper  with  phase  variation.  Later,  in  Section|^(beginning  on  page[4T])  on 
phase-only  tapers,  we  will  continue  to  report  @  but  with  the  average  amplitude 
gain  calculated  by  numerically  evaluating  the  integral  in  Q.  In  this  subcase  we 
have  no  Cauchy-Schwarz-Bunyakovsky  bound  on  the  values  this  average  mesa 
taper  loss  might  take. 

The  Airy  taper.  We  give  our  zero-phase  tapers  mesas  by  defining  pdf  p(f)  =  p(k^AB) 
to  be  constant  inside  a  circle  in  k  of  mesa  radius  kr  and  to  be  zero  outside  that  circle.  An 
inverse  2D  discrete-time  Fourier  transform  of  p(f)  and  subsequent  restriction  of  (pn  to  the 
array  aperture  then  gives  as  a  sampled  version 


=  2  jinc(27rA:r||BAn||) 


(7) 


of  the  Airy  pattern  from  optics,  using  a  jinc  function  defined  using  a  Bessel  function: 

for  X  ^  0, 


for  X  =  0. 


The  scaling  has  been  chosen  so  that  =  /square  7’(^)  di  =  1  as  required. 

Given  that  our  array  size  is  fixed,  kr  is  the  one  free  parameter  in  the  definition  of 
Average  mesa  taper  loss  is  presented  for  every  taper  in  this  report  for  which  a  mesa  is  an 
intended  design  feature,  and  for  all  amplitude  tapers  a  common  mesa  radius  parameter  kr  is 
used  in  calculating  Just  as  we  do  not  expect  a  filter’s  ripple  bandwidth,  3  dB  bandwidth, 
and  noise  bandwidth  to  be  the  same,  here  we  do  not  assume  that  kr  should  exactly  match 
the  design  radius  used  to  specify  the  mesa  size  in  the  optimization  constraints.  Instead  we 
determine  this  common  value  of  kr  by  computing  the  average  mesa  taper  loss  of  [Design 
Z4|for  a  range  of  values.  The  result  is  shown  in  Figure]^  The  lowest  average  mesa  taper 
loss  is  obtained  with  kr  =  0.0882/A,  so  that  value  is  used  for  calculating  average  mesa 
taper  loss  for  each  amplitude-taper  design  that  has  a  mesa. 

Of  course  taper  loss  is  indifferent  to  scaling,  so  rather  than  take  the  scaling  implied 
hy  h  =  ^,  we  can  if  we  wish  set  the  boresight  gain  to  approximately  unity.  The  scaling 
of  Fourier  pair  (pn  ^  p(f)  yields  ^  scale  (pn,  and  p(k^AB) 

by  |B^B|^/^7r(A:rA)^,  after  scaling  the  last  of  the  three  will  have  unit  height  in  the  mesa. 
After  scaling,  the  first  of  the  three  becomes  our  reference  taper  and  will  transform  to  an 
array  factor  of  approximately  unit  mesa  height  but  featuring  truncation-related  ripples  and 
sidelobes. 
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Figure  5:  Average  mesa  ta¬ 
per  loss  of  [Design  2!4|  versus 
tjjr,  the  angular  radius  of  the 
mesa.  We  choose  the  direction- 
cosine-space  mesa  radius  kr  = 
sin('^j.)/A  in  0  to  minimize 
this  [Design  2!4|  average  mesa 
taper  loss. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency)  13.80  dB 
max- weight  taper  loss  (max- weight  taper  efficiency)  29.09  dB 
average  mesa  taper  loss  (average  mesa  taper  efficiency)  0  dB 


(4.17%) 

(0.12%) 

(100%) 


Figure  6:  Reference  design  for  all  amplitude  tapers  with  mesas 
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Design  2.3  Weight-energy  minimization,  Taylor-like  pencil  beam 


The  first  set  of  optimizations  minimize  the  array-factor  energy  (or  weight  energy,  thanks  to 
Parseval).  Note  that  for  the  examples  with  mesas  this  is  not  exactly  the  same  as  minimizing 
the  weight-energy  taper  loss  since  the  boresight  gain  is  not  explicitly  fixed.  Nor  is  it  equiva¬ 
lent  to  minimizing  average  mesa  taper  loss,  because  the  average  mesa  gain  is  not  explicitly 
fixed.  This  can  lead  to  some  slightly  unexpected  results  when  comparing  designs,  as  shown 
later  in  the  discussions  for|Design  2.8|  and  [Design  2.12[  However,  the  peak  constraints  on 
the  mesa  do  approximately  fix  both  the  boresight  and  average  mesa  gains,  and  thus  the  net 
result  is  to  approximately  minimize  the  two  taper  losses  that  depend  on  total  weight  energy. 

This  design  is  a  pencil  beam  with  a  Taylor-like  response.  The  following  constraints  are 
applied  to  the  filter  (•  on  the  array-factor  plots): 


•  upperbound  the  array-factor  sidelobes  by  —35  dB  over  the  shelf  region,  from  2.5° 
outward. 

•  constrain  the  array  factor  at  boresight  to  0  dB. 

The  element  weights  are  nonnegative. 

The  points  where  constraints  are  applied  are  visible  on  each  of  the  array-factor  plots  (•). 
Due  to  the  symmetries  in  the  problem  formulation  (Figurej^,  these  points,  which  are  spread 
across  roughly  1/12  of  the  array-factor  period,  constrain  nearly  the  entire  universe.  Of 
course  if  we  were  to  remove  those  required  array-factor  symmetries  from  the  formulation, 
for  example  to  enable  the  inclusion  of  some  asymmetric  sidelobe  limits,  the  number  of 
variables  optimized  would  be  much  larger. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max- weight  taper  loss  (max- weight  taper  efficiency) 

weight-energy  level 


0.83  dB 
6.25  dB 
-35.71  dB 


(82.62%) 

(23.72%) 


Figure  7:  Weight-energy  minimization  of  Taylor-like  pencil  beam 
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Design  2.4  Weight-energy  minimization,  mesa  constraint 

Next,  we  test  the  mesa  constraint  on  our  design.  The  mainbeam  is  broadened  to  4.25°  so 
that  a  cluster  of  approximately  19  pencil  beams  will  fit  inside  the  mesa  (Figure]^.  The 
following  constraint  is  applied  to  the  filter  (•)  on  the  array-factor  plots): 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°. 


The  weight-energy  taper  loss  is  14.99  dB.  This  large  number  results  from  the  dissim¬ 
ilarity  between  the  present  wide  beam  and  the  reference  pencil  beam  of  [Design  2.1[  The 
average  mesa  taper  loss  is  a  more  appropriate  measure  because  it  compares  to  the  mesa- 
reference  wide  beam  of  [Design  2.2|  The  average  mesa  taper  loss  is  0.02  dB.  That  it  is 
very  small  is  no  surprise,  because  in  Figure]^  we  fixed  the  mesa  radius  of  [Design  2.2[to 
minimize  the  weight-energy  taper  loss  of  this  very  design.  The  element  weights  are  both 
negative  and  nonnegative. 

If  desired,  one  could  eliminate  the  hexagonal  rays  that  extend  from  the  mainbeam  by 
modifying  the  shape  of  the  circular  array. 


Figure  8:19  beams  packed  into 
a  hexagonal  cluster. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max- weight  taper  loss  (max- weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 
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14.99  dB 
29.97  dB 
0.02  dB 
-22.55  dB 


(3.17%) 

(0.10%) 

(99.47%) 


Figure  9:  Weight-energy  minimization  with  a  mesa  constraint 
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Design  2.5  Weight-energy  minimization,  mesa  constraint,  —35  dB  shelf 
constraint 

Here  we  combine  the  mesa  constraint  with  the  Taylor- like  shelf  constraint  to  get  a  broader 
beam  with  a  Taylor-like  response.  The  following  constraints  are  applied  to  the  filter  (•)  on 
the  array-factor  plots): 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°. 

•  upperbound  the  array-factor  sidelobes  by  —35  dB  over  the  shelf  region,  from  6.8° 
outward. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max- weight  taper  loss  (max- weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 
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(98.60%) 


Figure  10:  Weight-energy  minimization  with  a  mesa  constraint  and  a  —35  dB  Taylor-like 
shelf  constraint 
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Design  2.6  Weight-energy  minimization,  mesa  constraint,  —23  dB  shelf 
constraint 

Here,  we  change  the  shelf  height  from  —  35  dB,  a  reasonable  shelf  height  for  a  Taylor-like 
array,  to  —23  dB,  the  target  height  of  the  brim  for  our  sombrero  designs.  The  following 
constraints  are  applied  to  the  filter  (•)  on  the  array-factor  plots): 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°. 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  6.8° 
outward. 
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Figure  11:  Weight-energy  minimization  with  a  mesa  constraint  and  a  —23  dB  Taylor-like 
shelf  constraint 
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Design  2.7  Weight-energy  minimization,  Taylor-like  pencil  beam,  30° 
brim  constraints 

Here  we  go  back  to  a  pencil  beam  design  (we  remove  the  mesa  constraints)  and  add  the 
brim  constraints  (the  pedestal  and  shelf  constraints  in  Figure  [T])  to  produce  a  sombrero¬ 
shaped  array  factor  with  a  pencil-beam  mesa.  The  following  constraints  are  applied  to  the 
filter  (•)  on  the  array-factor  plots): 

•  constrain  the  array  factor  at  boresight  to  0  dB. 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  2.5° 
outward  (•)  on  the  array-factor  plots). 

•  lowerbound  the  array  factor  by  —25  dB  in  the  pedestal  region,  from  6.8°  to  30°  (o  on 
the  array-factor  plots). 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 

weight-energy  level 


8.13  dB 
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-28.41  dB 


(15.37%) 

(0.02%) 


Figure  12:  Weight-energy  minimization  of  a  Taylor- like  pencil  beam  with  30°  brim  con¬ 
straints. 
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Design  2.8  Weight-energy  minimization,  sombrero  constraints,  30°  brim 


Next,  mesa  constraints  are  reintroduced  to  generate  our  desired  sombrero  pattern.  The 
constraints  for  the  sombrero  pattern  are: 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°  (•)  on  the  array-factor  plots). 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  6.8° 
outward  (•)  on  the  array-factor  plots). 


•  lowerbound  the  array  factor  by  —25  dB  in  the  pedestal  region,  from  6.8°  to  30°  (o  on 
the  array-factor  plots). 


The  constraints  of  this  sombrero  design  include  all  of  the  constraints  from  the  Taylor- 
like  Mesa  design  ([Design  2^  plus  one  additional  (pedestal)  constraint.  Therefore  the  set 
of  possible  (feasible)  optimization  outcomes  of  this  sombrero  design  is  a  subset  of  the  set 
of  possible  optimization  outcomes  of  the  Taylor-like  mesa  design.  One  would  expect 

quality  (best ((small  set)))  <  quality  (best ((large  set))) 

if  (small  set)  C  (large  set).  It  appears  superficially  that  whether  this  expectation  is  met 
depends  on  the  quality  measure  chosen.  The 

...  ,  .  flb.lldB  Taylor-like  Mesa  design, 

(weight-energy  taper  loss>,B  =  I  Sombrero  design. 


and 


(energy)dB  = 


—22.43  dB  Taylor-like  Mesa  design, 
— 21.66  dB  Sombrero  design. 


The  weight-energy  taper  loss  results  do  not  match  our  expectations,  while  the  weight- 
energy  levels  do  meet  our  expectations.  The  energy  is  actually  minimized  by  the  opti¬ 
mization,  so  the  results  there  are  not  surprising.  Relationship  ([T])  can  be  written  in  decibel 
form  as  weight-energy  taper  loss  to  energy  is 


(weight-energy  taper  loss)dB  = 

(ideal  gain)dB  +  (energy)  dB  —  20  log^g  |  (boresight  array  factor)  | .  (8) 

Ideal  gain  is  not  affected  by  the  optimization,  but  the  boresight-array-factor  term  can  be. 
Here 

^  ,,  f  l.OOdB  Taylor-like  Mesa  design, 

-20  log„  I  (boresight  array  factor)  I  =  |  q  35  ^  s„„,brero  design. 

If  we  use  this  to  adjust  the  taper  loss. 


(weight-energy  taper  loss)dB 

—  20  logj^o  I  (boresight  array  factor)  | 


14.00  dB  Taylor-like  Mesa  design, 
14.88  dB  Sombrero  design. 


which  behaves  as  expected. 

We  conclude  from  this  array-factor  design  that  it  is  possible  to  optimize  a  sombrero 
pattern  for  an  array. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

weight-energy  level 


15.86  dB 
40.55  dB 
0.75  dB 
-21.68  dB 


(2.60%) 

(0.01%) 

(84.23%) 


Figure  13:  Weight-energy  minimization  with  sombrero  constraints  (30°  brim) 
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Design  2.9  Weight-energy  minimization,  sombrero  constraints,  45°  brim 

Here  is  another  example  of  the  sombrero  pattern,  with  the  brim  extending  out  to  45°.  The 
constraints  for  the  sombrero  pattern  are: 


•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°  (•)  on  the  array-factor  plots). 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  6.8° 
outward  (•)  on  the  array-factor  plots). 

•  lowerbound  the  array  factor  by  —25  dB  in  the  pedestal  region,  from  6.8°  to  45°  (o  on 
the  array-factor  plots). 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

weight-energy  level 


16.32  dB 
44.77  dB 
1.20  dB 
-21.22  dB 


(2.33%) 

(0.003%) 

(75.86%) 


Figure  14:  Weight-energy  minimization  with  sombrero  constraints  (45°  brim) 
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Design  2.10  Max- weight  optimization,  Taylor-like  pencil  beam 

The  second  set  of  optimizations  minimize  the  peak  magnitude  of  the  element  weights. 
From  @  we  can  obtain  a  relationship  analogous  to  ([8])  but  for  the  max-weight  taper  loss: 


(max-weight  taper  loss)dB  = 

(ideal  gain)dB  +  10  log^g  (A^I/imaxP)  -  20  log^g |  (Foresight  array  factor)  | .  (9) 

This  design  is  a  pencil  beam  with  a  Taylor-like  response.  The  following  constraints  are 
applied  to  the  filter  (•)  on  the  array-factor  plots): 

•  upperbound  the  array-factor  sidelobes  by  —35  dB  over  the  shelf  region,  from  2.5° 
outward. 


•  constrain  the  array  factor  at  Foresight  to  0  dB. 

Comparing  to  the  weight-energy  optimization  with  identical  constraints  ([Design  23]), 
note  that  the  peak  optimization  buys  us  nothing  (there  is  no  visible  change  in  the  peak 
coefficient)  and  costs  us  weight  energy. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency)  2.16  dB  (60.75%) 

max-weight  taper  loss  (max-weight  taper  efficiency)  5.32  dB  (29.35%) 

max-weight  level  —67.75  dB 


Figure  15:  Max- weight  optimization  of  a  Taylor- like  pencil  beam 
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Design  2.11  Max- weight  optimization,  mesa  constraint 

Next,  we  test  the  mesa  constraint  on  our  design.  The  mesa  constraints  broaden  the  main- 
beam  to  4.25°  to  fit  approximately  19  pencil  beams  inside  the  mesa  (Figure]^.  The  follow¬ 
ing  constraint  is  applied  to  the  filter  (•)  on  the  array-factor  plots): 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  Foresight 
to  4.25°. 

The  average  mesa  loss  of  this  case  is  significantly  worse  than  the  complimentary  weight- 
energy  optimization  case  ([Design  2'!4|).  Although  the  max- weight  taper  loss  is  reduced,  we 
also  have  very  large  and  undesirable  sidelobes  surrounding  the  mesa. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

max- weight  level 


22.58  dB 
22.62  dB 
8.07  dB 
-51.46  dB 


(0.55%) 

(0.55%) 

(15.58%) 


Figure  16:  Max-weight  optimization  with  a  mesa  constraint 
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Design  2.12  Max-weight  optimization,  mesa  constraint,  —35  dB  shelf 
constraint 

Here  we  combine  the  mesa  constraint  with  the  Taylor- like  shelf  constraint  to  get  a  broader 
beam  with  a  Taylor-like  response.  The  following  constraints  are  applied  to  the  filter  (•)  on 
the  array-factor  plots): 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°. 

•  upperbound  the  array-factor  sidelobes  by  —35  dB  over  the  shelf  region,  from  6.8° 
outward. 


Comparing  this  design  to|Design  2.5[  we  find  the  surprising  result  that  both  the  weight- 
energy  and  max-weight  taper  losses  have  improved,  albeit  marginally.  This  simply  reflects 
that  here  the  array  factor  peaks  up  at  boresight,  while  in  [Design  23]  the  array  factor  is 
depressed  at  boresight. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

max- weight  level 


14.37  dB 
28.99  dB 
0.13  dB 
-44.14  dB 


(3.66%) 

(0.13%) 

(97.05%) 


Figure  17:  Max- weight  optimization  with  a  mesa  constraint  and  a  —  35  dB  Taylor-like  shelf 
constraint 
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Design  2.13  Max- weight  optimization,  mesa  constraint,  —  23  dB  shelf 

constraint 

Here,  we  change  the  shelf  height  from  —  35  dB,  a  reasonable  shelf  height  for  a  Taylor-like 
array,  to  —23  dB,  the  target  height  of  the  brim  for  our  sombrero  designs.  The  following 
constraints  are  applied  to  the  filter  (•)  on  the  array-factor  plots): 

•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°. 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  6.8° 
outward. 


Here  again  the  improvement  in  max- weight  taper  loss  over  [Design  2.6|is  minimal  and 
accompanied  by  a  significant  loss  in  the  average  mesa  taper  loss. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

max- weight  level 


16.20  dB 
28.01  dB 
1.26  dB 
-46.07  dB 


(2.40%) 

(0.16%) 

(74.76%) 


Figure  18:  Max-weight  optimization  with  a  mesa  constraint  and  a  —  23  dB  Taylor-like  shelf 
constraint 
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Design  2.14  Max-weight  optimization,  Taylor-like  pencil  beam,  30°  brim 
constraint 


Here  we  go  back  to  a  pencil-beam  design  (we  remove  the  mesa  constraints)  and  add  the 
brim  constraints  (the  pedestal  and  shelf  constraints  in  Figure  [T])  to  produce  a  sombrero¬ 
shaped  array  factor  with  a  pencil-beam  mesa.  The  following  constraints  are  applied  to  the 
filter  (•)  on  the  array-factor  plots): 

•  constrain  the  array  factor  at  boresight  to  0  dB. 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  2.5° 
outward  (•)  on  the  array-factor  plots). 

•  lowerbound  the  array  factor  by  —25  dB  in  the  pedestal  region,  from  6.8°  to  30°  (o  on 
the  array-factor  plots). 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency)  11.45  dB 
max-weight  taper  loss  (max-weight  taper  efficiency)  30.22  dB 

max-weight  level  —42.86  dB 


(7.17%) 

(0.10%) 


Figure  19:  Max- weight  optimization  of  Taylor-like  pencil  beam  with  30°  brim  constraints 
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Design  2.15  Max-weight  optimization,  sombrero  constraints,  30°  brim 

Next,  mesa  constraints  are  reintroduced  to  generate  our  desired  sombrero  pattern.  The 
constraints  for  the  sombrero  pattern  are: 


•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°  (•)  on  the  array-factor  plots). 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  6.8° 
outward  (•)  on  the  array-factor  plots). 

•  lowerbound  the  array  factor  by  —25  dB  in  the  pedestal  region,  from  6.8°  to  30°  (o  on 
the  array-factor  plots). 

The  improvement  in  max-weight  taper  loss  over  [Design  2.8|is  modest. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

max-weight  level 


16.21  dB 
35.12  dB 
1.96  dB 
-37.96  dB 


(2.39%) 

(0.03%) 

(63.70%) 


Figure  20:  Max-weight  optimization  with  sombrero  constraints  (30°  brim) 
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Design  2.16  Max-weight  optimization,  sombrero  constraints,  45°  brim 

Here  is  another  example  of  the  sombrero  pattern,  with  the  brim  extending  out  to  45°.  The 
constraints  for  the  sombrero  pattern  are: 


•  bound  the  array  factor  between  —  1  dB  and  0  dB  in  the  mesa  region,  from  boresight 
to  4.25°  (•)  on  the  array-factor  plots). 

•  upperbound  the  array-factor  sidelobes  by  —23  dB  over  the  shelf  region,  from  6.8° 
outward  (•)  on  the  array-factor  plots). 

•  lowerbound  the  array  factor  by  —25  dB  in  the  pedestal  region,  from  6.8°  to  45°  (o  on 
the  array-factor  plots). 


As  in  many  of  the  previous  designs,  compared  to  [Design  2.9|we  see  a  small  improve¬ 
ment  in  the  max-weight  taper  loss  and  similar  increases  in  the  weight-energy  and  average- 
mesa  taper  losses. 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max-weight  taper  loss  (max-weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 

max-weight  level 


16.46  dB 
38.65  dB 
2.20  dB 
-34.43  dB 


(2.26%) 

(0.01%) 

(60.29%) 


Figure  21:  Max- weight  optimization  with  sombrero  constraints  (45°  brim) 
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3  Phase-Only  Tapers 

This  section  is  based  on  |l4||  but  includes  both  more  examples  and  material  on  taper  losses 
that  is  not  present  there. 

Earlier  taper  designs  in  this  report  all  featured  real  element  weights.  Optimization  of 
real  weights  or  even  complex  weights  is  practical  even  for  large  problems  because  the 
design  constraints  are  mathematically  convex:  given  any  two  tapers  and  that  satisfy 
the  constraints,  any  taper  obtained  by  linearly  interpolating  between  them,  like  .7  «„  +  -  3  Un 
or  .lun  +  .9un  or  Aun  +  .6un,  satisfies  those  constraints  as  well.  For  such  a  convex 
optimization  problem,  any  local  optimum  is  also  the  global  optimum.  Further,  there  are 
efficient  algorithms  for  finding  that  global  optimum. 

However,  robust  and  reliable  optimization  of  a  phase-only  array  taper  is  practical  only  in 
the  most  limited  of  special  cases,  simply  because  a  constant-amplitude  constraint  \hn\  =  1 
itself  is  nonconvex.  Interpolation  between  two  unit-magnitude  complex  numbers  does  not 
in  general  yield  a  unit-magnitude  result;  a  straight  line  between  two  points  on  the  circum¬ 
ference  of  a  circle  does  not  remain  on  that  circumference.  A  large  array  will  have  several 
thousand  of  these  nonconvex  constraints,  one  per  element  except  when  elements  are  related 
by  symmetry.  For  large  arrays  this  typically  gives  the  resulting  optimization  problem  a  vast, 
even  ghastly  number  of  local  optima,  and  there  is  no  reliable  way  of  discovering  which  of 
them  is  the  global  optimum. 

Viable  phase-only  strategies  then  tend  to  either  use  a  heuristic  method  to  pick  a  rea¬ 
sonable  suboptimal  design  or  to  use  a  preliminary  design  chosen  in  that  way  to  initialize 
a  follow-on  optimal  search  for  a  nearby  local  optimum.  Examples  of  heuristic  design  ap¬ 
proaches  can  be  found  in  |l5l|6l|7l[8l|9l[T0]].  Local  optimization  approaches  are  presented  in 
Illll[l2l[13[ll[l5l[l6l[13l.  Nonlocal  optimization  methods,  such  as  simulated  annealing 
and  genetic  algorithms,  have  also  been  used  IfTSl.  although  they  do  not  guarantee  global 
solutions  either.  There  are  also  hybrid  approaches  such  as  lfT9l.  in  which  a  small  number 
of  phase  basis  functions  is  heuristically  chosen  and  the  corresponding  coefficients  are  then 
locally  optimized.  Whether  any  of  these  procedures  yields  a  useful  result  depends  very 
much  on  the  particulars. 

A  simple  FM-chirp  heuristic  approach.  One  heuristic  approach  commonly  used  for 
uniform  line  arrays  in  ID  borrows  from  classical  FM  waveform  design.  The  core  idea  is 
simply  to  hope  that  the  magnitude  of  a  chirp’s  Fourier  spectrum  at  a  given  frequency  / 
is  proportional  to  the  amount  of  time  the  instantaneous  frequency,  given  for  chirp 
by  /inst(f)  =  spends  near  /.  It’s  very  much  like  an  integral  change  of  variable 

/  =  /inst(f)  or  t  =  /j~t(/)  so  that  dt  =  df.  The  chirp  spends  time  dt  between 

instantaneous  frequencies  /  and  /  +  df,  so  we  might  hope  that  its  Fourier  spectrum  will 

I  rlf —  ^  I 

have  magnitude  at  least  roughly  proportional  to  |  \  ■  The  messy  part  mathematically 

is  generally  dealing  with  the  inverse-function  aspect  of  this,  but  there  are  classic  cases  in 
which  the  relationships  are  not  difficult  to  work  out,  and  exploring  them  brings  out  the 
limitations  of  the  approach. 
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For  example,  for  f  G  f  >  f)  the  FM  chirp  waveform 


w{t) 


gj^(/3t)(t/T)  for  |f|  <  2"/2, 

0  otherwise. 


has  instantaneous  frequency  f3t/T,  which  sweeps  at  a  constant  rate  across  interval  (— f ,  f ) 
as  t  increases.  One  might  hope  then  that  its  Fourier  magnitude  spectrum  is  roughly  constant 
on  that  interval  and  zero  elsewhere.  In  fact  it  is  approximately  so  when  time-bandwidth 
product  l3T  is  large  and  only  then.  In  that  case, 


\w(f)\ 


some  constant  for  |/|  <p/2, 
0  otherwise. 


Use  unsampled  2D  tapers  and  their  transforms  for  simplicity.  We  begin  the  extension 
to  2D  by  establishing  a  simplified  notation.  Relating  that  new  notation  to  our  existing 
notation  is  a  tad  tedious  but  is  along  the  lines  of  the  usual  development  of  sampling  in 
signals-and-system  texts.  Let  the  orthonormal  columns  of  3  x  2  matrix  A  point  in  the 
azimuth  and  elevation  directions  when  looking  out  along  the  array  Foresight  vector.  If  two- 
vector  X  contains  azimuth  and  elevation  coordinates,  we  can  express  a  position  in  the  array 
plane  as  Ax,  and  we  can  express  the  array-plane  component  k  of  wavenumber  k  using  a 
two- vector  k  of  azimuth-elevation  coordinates  as  k  =  Afc. 

For  simplicity  let’s  then  assume  the  taper  is  obtained  by  sampling  a  continuous  func¬ 
tion  g{x)  at  the  values  of  x  corresponding  to  element  positions.  The  sampling  process 
can  be  represented  in  the  x  and  k  domains  and  continuous-”time”  (really  space)  Fourier 
transformation  in  2D,  where  each  sum  is  over  all  two-vectors  of  integers,  as 

9{x)  G{k) 

X  * 

A2|B^B|5j^(5(a;  -  A'^ABn)  Y^S{k  -  {{A^B)-^in) 

n  m 

II 

A2|B^B|^p(A^ABn)  5{x  -  A^ABn) 

n  ^  ^  n 

hn  ^  ^ 

H{k'^A^XB)  =H{k^XB). 

The  transform  in  the  last  line  follows  by  direct  evaluation  of  the  transform  integral.  The 
underbraces  identify  the  actual  taper  (left)  and  array  factor  (right)  in  terms  of  the  2D 
discrete-time  (space)  Fourier  pair  hn  ^  H{i)  used  earlier  for  amplitude  tapers.  Here 
g{x)  must  have  units  of  1/ (length)^  to  make  hn  dimensionless,  and  the  actual  array  factor 
Ff(k^AB)  =  H{\<PXB)  is  then  related  to  the  G{k)  that  we  will  use  here  by  k  =  Ak  and 

/f(fcWAB)  =  '^G{k-  i(A^B)-^m). 

m 

We  actually  have  no  further  need  of  such  details  of  scaling  and  periodic  replication,  as 
we  derive  results  only  up  to  a  scale  factor  and  assume  always  that  the  sampling-related 
replication  amounts  to  little  more  than  periodic  extension. 
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Brick-wall  beams  from  phase  bowls.  Extending  the  ID  chirp  idea  to  a  rectangle  in 
array-plane  coordinates  is  now  straightforward.  Beamwidth  and  time-extent  matrices  (3 
and  T  are  diagonal  with  positive  elements  on  those  diagonals: 


9{x) 


1 

some  area  ^ 

0 


for  T  G  <S, 
otherwise. 


(10) 


Here  set  <S  is  a  1  x  1  square  centered  on  the  origin.  In  2D  the  instantaneous  (really 
“instanspaceous”  since  x  comprises  position  coordinates,  but  we  lack  such  a  word)  fre¬ 
quency  of  a  signal  is  just  ki„st  =  ^V6*(cc),  with  the  gradient  taken  with  respect 
to  coordinate  vector  x.  The  phase  of  signal  g{x)  of  (fTO))  has  the  parabolic  “bowl”  shape 
0{x)  =  7rx^f3T~^x  and  instantaneous  frequency  k-mst  =  f3T~^x,  and  the  latter  ranges 
across  those  fcinst  for  which  /3“^fcinst  G  S.  And  indeed,  when  matrix  f3T  has  large  diagonal 
elements. 


|G(J;)| 


some  constant  for  /3  ^  fc  G  «S, 
0  otherwise. 


This  continuous  phase-only  taper  has  a  rectangular  outline  and  gives  an  approximately 
rectangular  beam  in  spatial  frequency  k  (sine  space). 

For  an  outline  elliptical  in  vector  x  we  could  use 


P(®) 


1  ^inx^/ST  ^x 
some  area 

0 


for  x^T  ^  cc  <  I , 
otherwise, 


where  2x2  symmetric  time  and  beamwidth  matrices  T  and  f3  here  are  positive  definite. 
The  major  and  minor  diameters  of  the  taper  are  given  by  the  eigenvalues  of  T.  Again 
fcinst  =  (3T~^x,  and  in  the  ellipse  where  g{x)  is  nonzero  it  ranges  over  the  ellipse  in  ki^st 
for  which  <  1/4,  so 


|G'(fc)| 


some  constant 
0 


for  ^k  <  1, 
otherwise. 


resulting  in  an  elliptical  beam  with  major  and  minor  diameters  given  by  the  eigenvalues  of 
f3.  Other  support  regions  can  be  used,  resulting  in  beams  with  matching  geometries. 

Any  of  these  continuous-space  functions  can  be  sampled  spatially  to  obtain  element 
weights  for  an  array.  The  primary  problem  with  such  “FM  beams”  is  that  the  approxima¬ 
tion  error  is  a  function  of  the  (eigenvalues  of  the)  space-beamwidth  product  /3T,  which 
(eigenvalues)  must  (all)  be  large.  Even  the  largest  radar  arrays  do  not  have  enough  ele¬ 
ments  for  this  approximation  to  be  reasonable — it  takes  hundreds  of  thousand  elements  to 
get  down  to  a  few  dB  of  ripple — so  we  are  reduced  to  tweaking  parameters  when  (each 
eigenvalue  of)  /3T  is  rather  small  in  hopes  of  “lucking  into”  a  decent  array  factor. 


Design  examples  approximating  brick-wall  beams.  We  illustrate  in  the  design  exam¬ 
ples  to  follow  by  presenting  several  actual  array  factors.  Our  array  and  desired  beam  here 
are  both  circular,  so  the  elliptical  case  above  applies.  That  circularity  makes  both  /3  and 
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T“^  into  scaled  identity  matrices,  so  we  can  parameterize  their  product  as  al  =  /3T~^ 
using  a  single  scalar  parameter  a.  Each  of  our  first  group  of  phase-only  tapers  therefore 
takes  the  form 

gjTralla; 


Our  spatial-extent  matrix  T  is  fixed  by  the  array  geometry,  so  choosing  our  one  available 
parameter  a  corresponds  to  choosing  beamwidth  matrix  {i.  Given  a  desired  beam  size, 
we  have  no  degrees  of  design  freedom  left.  Here  we  simply  pick  several  representative  a 
values  for  the  sake  of  illustration. 

We  do  notice  that  the  boresight  value  of  the  array  factor  seems  especially  sensitive  to  a, 
so  given  a  general  range  for  a  that  gives  us  roughly  the  beamwidth  we  seek,  we  generally 
end  up  fine  tuning  a  to  set  this  boresight  value  to  something  nonextreme. 


for  1 1  a?  1 1  in  our  4507-element  array, 
otherwise, 


(11) 


A  crude  attempt  at  a  Gaussian  beam.  The  phase-bowl  design  examples  just  discussed 
have  instantaneous  frequencies  fcinst  that  are  linear  in  position  x  and  so  implicitly  aim  at 
having  their  array  factors  approximate  “brick  wall”  spectra.  It  is  no  surprise  therefore  that 
nasty  array-factor  ripples  arise  as  a  result,  and  it  certainly  seems  reasonable  to  attempt  some 
sort  of  gradual  rolloff  as  an  alternative,  in  hopes  that  the  approximations  will  be  better,  that 
ripple  magnitudes  will  be  significantly  lower. 

Towards  this  end,  consider  the  taper 


9{x) 


for  ||a;||  <  Tq, 
otherwise. 


(12) 


Different  authors  define  the  error  function  with  different  scalings,  but  here  z  =  erf 
refers  to  the  odd  function  defined  for  positive  arguments  by 

y  =  eri{z)^^[  e~^^du. 

Jo 

A  little  algebra  shows  instantaneous  frequency  as  a  function  of  x  to  be 

fcinst(®)  =  ^Vphase(a;)  =  ^^x  ^^^eTi~\\\x\\/ro), 
and  change  of  variable  k  —  ki„stix)  then  results  (shown  in  lIU)  in  the  ratio  of  differential 


areas 


Then 


^feinst(fe)  ^  M  ^erf(||fc||/v/^)\  ^_||fc||2/(2fc2) 
dk  tt/c^  V  Ilfcll/A/^  / 


k\\/^  j 


(13) 


|G(lfe)|  »  (some  constant)  x  ,14, 

at  least  for  sufficiently  high  roko,  which  here  serves  the  role  of  a  time-bandwidth  product. 
The  quantity  in  large  parenthesis  in  ([T3])  goes  to  unity  as  k  goes  to  zero,  and  the  error 
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function  approaches  i/2/7r  for  large  arguments,  so  the  Fourier  spectrum  approximated  in 
(fT4|)  is  essentially  Gaussian  for  ||fc||  <C  ko  but  rolls  off  somewhat  faster  than  Gaussian  for 
||fc||  >fco. 

We  present  designs  below  for  ro  =  20.31A,  and  several  different  values  of  k^.  It  is  no 
surprise  that  the  array  factor  approximates  the  value  here  more  closely  when  tq  k^  is  large 
and  that  for  beam  sizes  of  interest  r^k^  in  fact  is  fairly  small.  Still  the  resulting  array- 
factor  ripple  magnitudes  can  sometimes  be  made  far  lower  than  for  the  phase-bowl  tapers 
by  choosing  ro  and  fco  carefully.  Generally  approximation  quality  is  far  more  sensitive  to  fco 
than  to  ro,  and  again  the  value  of  the  array  factor  at  Foresight,  at  the  origin,  is  particularly 
sensitive.  Very  small  deviations  from  the  k^  values  shown  can  triple  or  quadruple  mainbeam 
ripple. 


Sombrero-shaped  array  factor.  One  possible  way  to  generate  a  sombrero- shaped  phase- 
only  array  taper  is  to  replace  ||a?||  in  (fTl))  with  some  function  /(||a;||)  giving 


9{x) 


for  ||a;||  <  ro, 
otherwise. 


(15) 


We  present  several  designs  below  using  /(||a;||)  =  6||a;||  -I-  (o"*  -I-  —  a.  As  with 

the  previous  phase-only  designs,  small  deviations  in  any  one  of  the  five  parameters  k^,  ro, 
a,  b,  or  m  cause  major  changes  in  the  array  factor. 


Presentation  of  individual  designs.  For  each  phase-only  design  that  follows,  the  array- 
factor  magnitudes  are  shown  normalized  to  a  maximum  magnitude  of  0  dB.  Array-factor 
phases  are  nonzero  but  are  not  shown.  The  phase  of  the  element  weights  are  plotted  in 
3D  in  the  lower  right.  The  plot  above  that  represents  instantaneous  frequency  along  an 
arbitrary  cut  through  the  array  origin  with  fcinst  equal  to  the  plotted  quantity  times  a  unit 
vector  in  the  cut  direction.  The  average  mesa  taper  loss  of  @  (page[T^  is  calculated  using 
^  of  0  (page  with  a  mesa  radius  parameter  kr  customized  to  that  specific  design  by 
searching  for  a  minimum  in  the  spirit  of  Figure  |5]  (page  [TT]). 
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Design  3.1  Linear-FM  chirp,  small  a 

Here  we  use  LFM  chirp  ([TT])  with  a  small  a.  The  3  dB  beamwidth  of  this  design  is  similar  to 
the  beamwidth  of  our  amplitude-taper  pencil-beam  designs  ([Design  23]  and  [Design  2.10|), 
but  the  rolloff  in  this  design  is  much  slower. 
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Results: 
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average  mesa  taper  loss  (average  mesa  taper  efficiency)  1 . 1 1  dB  (77.36%) 

Figure  22:  Linear-FM  chirp,  small  a 
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Design  3.2  Linear-FM  chirp,  medium  a 

Here  we  use  LFM  chirp  ([TT])  with  a  medium  a.  The  3  dB  beamwidth  of  this  design  is  similar 
to  the  beamwidth  of  our  amplitude-taper  mesa  constrained  designs  ([Design  2!4|  and  [Design 
[2.11[).  Compared  to  those  mesa  constrained  designs,  the  ripple  in  this  design  is  predictably 
bad,  the  rolloff  is  much  slower  and  the  overall  sidelobe  level  is  higher.  However,  we  do  see 
a  significant  improvement  in  the  max-weight  taper  loss  when  compared  to  our  optimized 
amplitude-taper  designs. 
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Figure  23:  Linear-FM  chirp,  medium  a 
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Design  3.3  Linear-FM  chirp,  large  a 

Here  we  use  LFM  chirp  ([TT])  with  a  large  a,  which  gives  us  a  design  with  a  larger  mesa. 
Both  the  ripple  in  the  mesa  and  the  difference  in  the  boresight  array  factor  are  very  sensitive 
to  a. 
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Figure  24:  Linear-FM  chirp,  large  a 
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Design  3.4  Nonlinear-FM  chirp  for  narrow  Gaussian  beam 

Here  we  use  (fTl])  with  a  fixed  tq  and  a  small  k^. 
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Results: 

weight-energy  taper  loss  (weight-energy  taper  efficiency)  1 1 . 15  dB  (7.68%) 
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average  mesa  taper  loss  (average  mesa  taper  efficiency)  1.25  dB  (74.98%) 

Figure  25:  Nonlinear-FM  chirp  for  Gaussian  beam,  small 
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Design  3.5  Nonlinear-FM  chirp  for  medium- width  Gaussian  beam 

Here  we  use  (fTl])  with  a  fixed  tq  and  a  medium  k^.  This  design  is  comparable  to  [Design 
|3.2|  and  our  amplitude-taper  mesa  designs  ([Design  2'!4|  and  [Design  iTH).  Compared  to 
Design  2.^the  weight-energy  and  average-mesa  taper  losses  are  both  somewhat  worse,  and 
as  expected  the  max- weight  taper  loss  is  significantly  improved.  Compared  to  [Design  2.11 
all  the  taper  losses  are  significantly  improved.  The  ripple  in  the  mesa  is  much  smaller  in 
this  design  than  in  [Design  3.2[and  the  roll-off  and  sidelobe  levels  are  similar. 
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Figure  26:  Nonlinear-FM  chirp  for  Gaussian  beam,  medium  ko 


55 


Design  3.6  Nonlinear-FM  chirp  for  wide  Gaussian  beam 

Here  we  use  (fTl])  with  a  fixed  tq  and  a  large  Uq.  These  parameters  generate  a  design  with 
a  larger  mesa  comparable  to  [Design  3!^  The  ripple  in  the  mesa  is  much  smaller  in  this 
design  than  in  [Design  3.3[  but  the  roll-off  and  sidelobe  levels  are  not  as  good. 
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Results: 

weight-energy  taper  loss  (weight-energy  taper  efficiency)  24.94  dB  (0.32%) 

max-weight  taper  loss  (max-weight  taper  efficiency)  24.94  dB  (0.32%) 

average  mesa  taper  loss  (average  mesa  taper  efficiency)  0.84  dB  (82.46%) 

Figure  27:  Nonlinear-FM  chirp  for  Gaussian  beam,  large  /cq 
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Design  3.7  Nonlinear-FM  chirp  for  sombrero  beam,  first  try 

Here  we  use  For  this  design,  the  mesa  is  flat  for  approximately  2.5°,  there  is  a  — 15  dB 

shelf  to  approximately  30°  and  a  sidelobe  level  of  —30  dB  at  90°. 
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Figure  28:  Nonlinear-FM  chirp  for  sombrero  beam,  first  try 
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Design  3.8  Nonlinear-FM  chirp  for  sombrero  beam,  better  transition 

Here  is  another  design  using  (fTS)). 
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Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max- weight  taper  loss  (max- weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 


15.32  dB  (2.94%) 
15.32  dB  (2.94%) 
1.37  dB  (72.95%) 


Figure  29:  Nonlinear-FM  chirp  for  sombrero  beam,  better  transition 
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Design  3.9  Nonlinear-FM  chirp  for  sombrero  beam,  pointy  beam 

Here  we  use  This  design  has  a  pointy  beamed  mesa  with  a  slow  transition  to  a  shelf 
and  rolloff  comparable  to  [Design  3.7 


62 


OdB 
-10  dB 
-20  dB 
-30  dB 
-40  dB 
-50  dB 
-60  dB 


0-3  dB 
0-2  dB 

IdB 

OdB 

-IdB 

-2dB 

-3dB 


galaxy 

zoom 


Parameters: 


world 

zoom 


natiomi 

zoom 

mdB 


nation 

zoom 

bdB 


dB  slice  plot 


instantaneous 
frequency 
vs.  r 


Phase 

weights 


A:o  =  0.3291/A 
ro  =  24A 
a=  12A 
5  =  0.68 
m  —  Q 

kr  =  0.5569/A 


Results: 


weight-energy  taper  loss  (weight-energy  taper  efficiency) 
max- weight  taper  loss  (max- weight  taper  efficiency) 
average  mesa  taper  loss  (average  mesa  taper  efficiency) 


15.78  dB  (2.64%) 
15.78  dB  (2.64%) 
1.35  dB  (73.36%) 


Figure  30:  Nonlinear-FM  chirp  for  sombrero  beam,  pointy  beam 
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4  Conclusions 


optimizing  sombrero- shaped  amplitude  (or  complex)  tapers  is,  although  not  a  trivial  prob¬ 
lem,  one  that  we  can  successfully  solve  with  our  current  optimization  tools.  Of  the  two 
amplitude-taper  optimization  objectives  presented,  minimizing  total  transmitter  power  while 
(approximately)  fixing  boresight  or  mesa-averaged  radiated  power  density  generally  pro¬ 
duces  better  results  than  minimizing  the  maximum  element  power.  For  the  most  part  the 
latter  simply  increases  the  weight-energy  and  average  mesa  taper  losses  without  bringing 
the  max-weight  taper  loss  anywhere  near  a  reasonable  number  and  at  the  additional  cost 
of  irregular  and  sometimes  large  sidelobes.  There  may  be  specific  situations  where  some 
gentle  constraints  on  the  peak  element  weights  would  be  useful,  but  in  general  if  the  array 
to  be  designed  is  heavily  peak  limited  then  an  amplitude  taper  appears  to  be  a  rather  poor 
choice.  Unfortunately,  current  digital-array  driver  concepts  are  inherently  peak  limited. 

Phase-only  tapers  do,  when  designs  comparable  to  optimized  amplitude  tapers  are  dis¬ 
covered,  significantly  improve  the  taper  loss  in  the  peak-limited  element  case.  However,  we 
cannot  generate  a  sombrero  shaped  phase-only  design  with  the  same  degree  of  control  and 
performance  as  an  amplitude  taper  when  the  same  array  geometry  is  used.  The  phase-only 
sombrero  shaped  designs  shown  here  were  generated  by  tweaking  five  different  parameters 
in  (fTS)).  Each  of  these  five  parameters  has  significant  effects  on  the  array  factor,  but  those 
effects  are  not  systematic.  The  designs  were  a  result  of  guessing  a  starting  point  and  hand 
adjusting  the  parameters  to  attempt  to  design  a  reasonable  sombrero  shaped  array  factor. 
This  process  could  be  formalized  by  directly  optimizing  such  parameters  or  even  the  phases 
themselves,  but  the  presence  of  large  numbers  of  local  minima  would  make  it  impractically 
expensive  to  exhaustively  search  for  the  global  minimum. 
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Design 

Optimization 

variable  level 

weight- 

energy 

taper 

loss 

max- 

weight 

taper 

loss 

average 

mesa 

taper 

loss 

Uniform 

N/A 

N/A 

OdB 

OdB 

N/A 

Taylor  Pencil 

weight  energy 

-35.71  dB 

0.83  dB 

6.25  dB 

N/A 

Mesa 

weight  energy 

-22.55  dB 

14.99  dB 

29.97  dB 

0.02  dB 

Taylor  Mesa  —35  dB 

weight  energy 

-22.43  dB 

15.11  dB 

30.98  dB 

0.06  dB 

Taylor  Mesa  —23  dB 

weight  energy 

-22.54  dB 

15.00  dB 

30.20  dB 

0.02  dB 

Taylor  Brim 

weight  energy 

-28.41  dB 

8.13  dB 

37.01  dB 

N/A 

Sombrero  30° 

weight  energy 

-21.68  dB 

15.86  dB 

40.55  dB 

0.75  dB 

Sombrero  45° 

weight  energy 

-21.22  dB 

16.32  dB 

44.77  dB 

1.20  dB 

Taylor  Pencil 

max  weight 

-67.75  dB 

2.16  dB 

5.32  dB 

N/A 

Mesa 

max  weight 

-51.46  dB 

22.58  dB 

22.62  dB 

8.07  dB 

Taylor  Mesa  —35  dB 

max  weight 

-44.14  dB 

14.37  dB 

28.99  dB 

0.13  dB 

Taylor  Mesa  —23  dB 

max  weight 

-46.07  dB 

16.20  dB 

28.01  dB 

1.26  dB 

Taylor  Brim 

max  weight 

-42.86  dB 

11.45  dB 

30.22  dB 

N/A 

Sombrero  30° 

max  weight 

-37.96  dB 

16.21  dB 

35.12  dB 

1.96  dB 

Sombrero  45° 

max  weight 

-34.43  dB 

16.46  dB 

38.65  dB 

2.20  dB 

Table  1 :  Amplitude-taper  results 
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weight- 

max- 

average 

energy 

weight 

mesa 

taper 

taper 

taper 

Design 

variable 

loss 

loss 

loss 

Linear-FM  chirp 

=  0.003/A2 

6.45  dB 

6.45  dB 

1.11  dB 

Linear-FM  chirp 

a  =  0.006/A2 

15.04  dB 

15.04  dB 

0.84  dB 

Linear-FM  chirp 

a  =  0.015/A2 

29.46  dB 

29.46  dB 

0.58  dB 

Nonlinear-FM  chirp 
Gaussian  beam 

ko  =  0.07/A 

11.15  dB 

11.15  dB 

1.25  dB 

Nonlinear-FM  chirp 
Gaussian  beam 

fco  =  0.1068/A 

16.04  dB 

16.04  dB 

1.07  dB 

Nonlinear-FM  chirp 
Gaussian  beam 

fco  =  0.3291/A 

24.94  dB 

24.94  dB 

0.84  dB 

fco  =  0.3291/A 

Nonlinear-FM  chirp 
Sombrero  beam 

To  =  24A 
a  =  12.048A 
h  =  0.697 
m  =  5 

18.26  dB 

18.26  dB 

1.32  dB 

fco  =  0.3291/A 

Nonlinear-FM  chirp 
Sombrero  beam 

ro  =  20.31A 
a  =  lOA 

6  =  0.4 
m  =  6 

15.32  dB 

15.32  dB 

1.37  dB 

fco  =  0.3291/A 

Nonlinear-FM  chirp 
Sombrero  beam 

ro  =  24A 
a  =  12A 
h  =  0.68 

15.78  dB 

15.78  dB 

1.35  dB 

m  =  6 


Table  2:  Phase-only  taper  results 
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